Quasiparticle vanishing driven by geometrical frustration 



O 

o 

(N 



i 

-4-* 



i 

c 

o 
o 



> 

O 
m 
O 

03 



O 

o 



X 



A. E. Trumper, C. J. Gazza, and L. 0. Manuel 
Instituto de Ftsica Rosario ( CONICET) and Universidad Nacional de Rosario, 
Boulevard 27 de Febrero 210 bis, (2000) Rosario, Argentina 
(Dated: February 2, 2008) 

We investigate the single hole dynamics in the triangular t — J model. We study the structure 
of the hole spectral function, assuming the existence of a 120° magnetic Neel order. Within the 
self-consistent Born approximation (SCBA) there is a strong momentum and t sign dependence of 
the spectra, related to the underlying magnetic structure and the particle-hole asymmetry of the 
model. For positive t, and in the strong coupling regime, we find that the low energy quasiparticle 
excitations vanish outside the neighbourhood of the magnetic Goldstone modes; while for negative 
t the quasiparticle excitations are always well defined. In the latter, we also find resonances of 
magnetic origin whose energies scale as (J/t) 2 ^ 3 and can be identified with string excitations. We 
argue that this complex structure of the spectra is due to the subtle interplay between magnon- 
assisted and free hopping mechanisms. Our predictions are supported by an excellent agreement 
between the SCBA and the exact results on finite size clusters. We conclude that the conventional 
quasiparticle picture can be broken by the effect of geometrical magnetic frustration. 



I. INTRODUCTION 

The study of doped antiferromagnets has become a 
central topic in the area of strongly correlated electrons 
systems, mainly after Anderson's proposition^ about the 
probable existence of non-Fermi liquid behavior once a 
Mott insulating state is hole-doped. A valuable exper- 
iment that gives access to single-particle properties of 
many-electron systems is the angle resolved photoemis- 
sion spectroscopy (ARPES). Even though in the last 
years there has been a rapid pace of developments in 
ARPES experiments, there are still controversies in many 
strongly interacting cases, like the cuprates, where the 
difficulty to distinguish a coherent <5-function peak from 
an incoherent part of the spectra does not allow a proper 
interpretation of the available data 2 *^. A useful mi- 
croscopic model widely used to describe the low-energy 
physics involved in photoemission experiments of anti- 
ferromagnetic Mott insulators^ is the t — J model on the 
square lattice^. 

Attention has also been concentrated on triangular 
Mott insulators. These systems are attractive due to 
the presence of both, strong correlation and geometri- 
cal magnetic frustration, leading to non-conventional be- 
havior. Experimental realizations of such systems are 
Cs^CuCl/^, the k family of the organic charge trans- 
fer salts (BEDT - TTF) 2 X 7 , and the silicon surfaces 
K/Si (111):B and SiC(OOOl) 8 ' 9 . In a recent publication, 
Koretsune and Ogata— have studied the effect of finite 
doping in the triangular Heisenberg model. They found 
that a resonating-valence-bond (RVB) ground state is fa- 
vored when t is positive, while a Nagaoka's ferromagnetic 
ground state is stabilized for negative t. More recently, 
the role of the t sign and frustration has been addressed 
in the superconductor Na x Co02 where the Co atoms 
form the triangular latticed. These remarkable issues 
lead us to investigate the spectral function of a hole in 
a triangular Mott insulator, so as to analize the conse- 
quences of the geometrical frustration and the t sign on 



the hole motion. In particular, our aim is to study the 
existence of quasiparticle (QP) excitations and higher- 
energy features in the spectral functions. To this pur- 
pose, we focus on the hole dynamics in the triangular 
t—J model, and we solve it using the Lanczos exact diag- 
onalization method and the self-consistent Born approx- 
imation (SCBA). In the non- frustrated case, it is known 
that the SCBA reproduces quite well the exact diagonal- 
ization results on small clusters^ and the angle-resolved 
photoemission spectroscopy (ARPES) experiments^. In 
particular, the SCBA predicts the existence of a well de- 
fined quasiparticle peak (magnetic polaron) at low ener- 
gies for the whole Brillouin zone (BZ). Furthermore, the 
SCBA is able to capture finite lifetime resonances above 
the QP peaks. These resonances have been explained by 
the string picture since their energies scale as (J/t) 2 / 3 , 
as it happens in the Ising limit. In our triangular frus- 
trated case, the presence of the 120° Neel order gener- 
ates two mechanisms for hole motion in the SCBA effec- 
tive Hamiltonian, one being tight-binding like (or free) 
and the other, magnon-assisted. We will show that the 
t sign tunes the interference between both mechanisms, 
producing pronounced differences in the low energy as 
well as in the higher energy region of the spectral func- 
tions. For the strong coupling regime, we find an excel- 
lent agreement between the SCBA and the exact results 
on a 21 sites cluster. In the thermodynamic limit, the 
SCBA spectral functions present a rich structure. No- 
tably, for t < a quasiparticle excitation is always de- 
fined for all k of the Brillouin zone, while for t>0 in the 
strong coupling regime the quasiparticle signal survives 
only around the Goldstone magnetic wave vectors. By 
computing the quasiparticle energy scaling with J / 1 it is 
possible to get some insight of the magnetic cloud sur- 
rounding the holei 3 ^. We argue that for t < the local 
environment around the hole enhances its AF character, 
while for t > the ferromagnetic character is favored. 
We find that the single-hole dynamics in the thermody- 
namic limit behaves in the opposite way as in very small 
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clusters, like the three-sites problemiS. 

Furthermore, for small values of J/t and t < 0, higher- 
energy resonances above the quasiparticle peak appear 
showing a (J/t) 2 / 3 energy scaling. Although the triangu- 
lar t — J z model is completely different to the square lat- 
tice case, we continue identifying these resonances with 
string- like excitations. For t>0, instead, we do not ob- 
serve higher energy resonances that can be identified with 
strings. This behavior observed at higher energies can be 
also traced back to the particular features of the magnetic 
environment around the hole mentioned above. 

The paper is organized as follows. In section II we 
obtain the effective Hamiltonian and the one hole prop- 
erties within the SCBA. In section III we present the 
results: (HIA) the comparison between exact and SCBA 
results in finite size clusters, (HIB) the quasiparticle fea- 
tures, and (IIIC) the general description of the spectral 
functions. Finally, in section IV we draw the concluding 
remarks. 



II. T-J MODEL AND SELF CONSISTENT 
BORN APPROXIMATION 



The t— J model is defined as 



H = -*E(Ca* + h - c -) + J E s < ■■ • ^ . (!) 



(id) 



(id) 



where the first term represents the hopping between near- 
est neighbors (< i,j >) of the triangular lattice, with the 
constraint of no double occupancy, Cj i(T = Ci )(T (l — n^-a-), 
and the second term represents the antifcrromagnetic 
(AF) Heisenberg part with J>0. 

The study of the single hole motion in an AF requires 
a correct description of the magnetic low energy sector 
since the injected hole will couple strongly with these spin 
excitations. In the literature there is firm numerical 14 
and analytical^ evidence that the triangular Heisenberg 
antiferromagnet has a 120° Neel order ground state, and 
that spin wave theory reproduces quite well its low energy 
behavior' 1 - 6 . So that, we assume an ordered ground state 
with a magnetic wave vector Q = (4tt/3, Q) lying in the 
x—z plane, and spin waves as the low energy excitations. 
In order to simplify the calculation we perform a unitary 
transformation to local spin quantization axis -primed 
operators- in order to have a ferromagnetic ground state 
in the z' direction. Then, we use the spinless fcrmion 



representation for electrons c- T = h\, 



At' _ 



hiS, , and 



the Holstein-Primakov representation for spin operators 
(to order 1/3) Sf' ~ \{a\ + a^Sf' ~ |(aj - a t ),Sf = 
h — a\a,i, where h\ creates a spinless hole at site i, and 
a\ creates a spin deviation. 

We replace the above expresions in eq. (JTJ and, af- 
ter Fourier and Bogoliubov transformations, the effective 
Hamiltonian results 
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with e k = — i7k and w q = §J>/(1 - 37 q )(l + 67 q ), the 
bare hole and magnon dispersions, respectively. 



Mcq = «(/?kU-q - /?k-qWq) 

is the bare hole-magnon vertex interaction with the geo- 
metric factors 7k = | J2 e cos(k.e) and /3k = Y^e sm (k.e) 
(e's are the positive vectors to nearest neighbors), and 
u q and u q are the usual Boguliubov coefficients. 
In the Hamiltonian J2J -constant terms related to mag- 
netic energy are omitted- the free hopping hole term 
implies a finite probability of the hole to move without 
emission or absorption of magnons. This is a direct con- 
sequence of the underlying non-collinear magnetic struc- 
ture. The hole-magnon interaction adds another mech- 
anism for the charge carrier motion which is magnon- 
assisted. The latter is responsible for the spin-polaron 
formation when a hole is injected in a non-frustrated 
antiferromagnetiSiii. We will show the existence of a 
subtle interference between both processes that turns out 
to be dependent on the t sign. An important quantity 
that allows us to study the interplay between such pro- 
cesses is the retarded hole Green function that is defined 
as G£(w) = {AF\h k (u+i ^_ H) hl\AF), where \AF) is the 
undoped antiferromagnetic ground state with a 120° Neel 
order. In the SCBA \AF) is the spin wave AF ground 
state, and after neglection of the crossing diagrams, the 
Dyson's equation for the self energy at zero temperature 
results 
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We have solved numerically this self-consistent equa- 
tion for Sk(w), and calculated the hole spectral func- 
tion Ak(w) = — jlmGj(w), and the quasiparticle (QP) 

weight Zk = (l — ^%t^) \E k , where the QP energy is 
given by E k = S k (£ k ). 



III. RESULTS 

A. Exact and SCBA comparison 

We have computed the one hole spectral function 
Ak(w) using the Lanczos exact diagonalization method 
(ED) and the SCBA on finite size clusters. To test the 
reliability of the SCBA predictions we have compared 
both techniques on the N — 21 sites cluster, which pre- 
serves the spatial symmetries of the infinite triangular 
lattice^. In general, we have found a good agreement 
between ED and SCBA results for all momenta of the 
BZ and for J/\t\ < 1. In Fig. ^we show the spectral 
functions at three representative k points of the BZ (see 
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FIG. 1: Spectral functions versus frequency for J/\t\ = 0.4 
and N = 21 for momenta T = (0,0), A = ^(-1,3V3) and 
B — ^(2,V3)- k shown as filled circles in the inset (the 
crosses represent the other momenta). The upper (lower) 
panel displays the results for positive (negative) t. The solid 
and dashed lines are the exact and SCBA results, respectively. 
rj = 0.1 and ui = 10000 has been used. 



inset Fig. ^| in the strong coupling regime, J/\t\ = 0.4, 
and for both signs of the transfer integral t. 

The exact structure of the spectral functions strongly 
depends on the momentum because the hole strongly 
scatter off antiferromagnetic fluctuations^. Furthermore, 
the geometrical frustration is manifested in a marked de- 
pendence of the spectral function with the t sign. Re- 
markably, the SCBA reproduces very well these exact 
results in the whole range of energies. In some cases 
from the scale of the figure it cannot be clearly distin- 
guished a quasiparticle peak at low energies. By decreas- 
ing r] + appropriately we have checked, however, the ex- 
istence of a QP peak for all k of the BZ. In particular, 
for negative (positive) t the k = A (k = T) point be- 
comes the momentum of the quasiparticle ground state. 
The energy of the hole can be splitted in two contri- 
butions: E h = E^-i - E% s = E? oc + E h deloc , where 
E'loc — ^*N-i ~ 1S the energy of the localized hole, 
and E% eloc = Em-i — E t ^ 1 is the energy gain due to 
delocalization of the hole— As the SCBA does not take 
into account the localized hole energy we have shifted the 
exact energy spectra by an amount E^ oc = —0.5686 J for 
A = 21 in order to make a proper comparison with the 
SCBA results. The agreement between both techniques 
is similar to that found in the 4x4 square lattice besides 
the presence of stronger antiferromagnetic fluctuations in 
our frustrated case. It was shown that this agreement is 
a consequence of the small leading nonzero contribution 
to the self energy from the crossing diagrams 19 . 



B. Quasiparticle 

The momentum and the t sign dependence of the 
spectra aforementioned is strengthened in the thermo- 
dynamic limit. In particular, we observe important dif- 



ferences in the quasiparticle weight z^ through the BZ 
under the change of t sign. Using lattice sizes up to 
N = 2700 we have extrapolated the results to the ther- 
modynamic limit. In Fig|2|we show the values of z\^ as a 
function of J/t for representative points along high sym- 
metry axes of the BZ . The most salient feature we have 
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FIG. 2: Quasiparticle weight vs J/t. Upper (lower) panel is 
for t > (t < 0). The location of the F, M' , K, M is displayed 
in Fig. El 

obtained is the vanishing of the quasiparticle weight at 
some momenta for positive t. At the M and M' points 
z w — below J/t ~ 2.5 and 1.5, respectively. It is worth 
noticing the robustness and then the rapid decay to zero 
of zm 1 as J/t — > 1.5. A similar behavior has the QP sig- 
nal at r, the ground state momentum, but it goes to zero 
only when J/t — > 0. Finally, the QP signal is very weak, 
but finite, at the magnetic wave vector K for all the J/t 
range studied; for instance, z^ ~ 0.008 when J/t = 10. 
Surprisingly, for J >> t, zk does not seem to approach 
one as it turns out in the square lattice casein . The non 
existence of quasi-particle excitations is a striking mani- 
festation of the strong interference between the free and 
magnon-assisted hopping processes, tuned by the t sign. 
For the strong coupling regime and a large region of the 
BZ, the coherence of the QP is lost because of the quite 
different group velocities of the spin fluctuation and the 
tight binding hole motion, suggesting the likelihood of 
a spin-charge separation scenario 1,3 '— for the triangular 
t — J model. 

On the other hand, for negative t the quasiparticle 
weight is finite for all momenta and for J/t ^ 0. This 
behavior is similar to the one found in the non-frustrated 
case— It is interesting to note that now the QP ground 
state momentum is M (K) for J/\t\ < 1.2 (J/\t\ > 1.2). 
The spectral weight of the QP becomes the most robust 



FIG. 3: Triangular Brillouin zone. The shaded areas cor- 
respond to the regions where the QP weight is finite for 
J/\t\ = 0.4 and positive t. 



for the ground state momenta. The interchange in the 
k dependence of the QP weight, evidenced in a compar- 
ison between both panels of Fig. [3 could be thought as 
a remnant of the particle-hole symmetry that shifts the 
momenta by Q under t sign reversal. 

In Fig. the shaded areas correspond to the regions 
where the QP weight is finite for J/\t\ = 0.4 and positive 
t. There is no QP signal for momenta outside the neigh- 
borhood of the magnetic Goldstone modes (k = 0,Q). 
The presence of QP signal in the non-equivalent ver- 
tices of the BZ is an artifact of the used approxima- 
tion: the spin wave magnetic ground state has a defi- 
nite chirality (Q), but our vertex interaction does not 
discriminate between ±Q, resulting the whole spectral 
function degenerate in both momenta. The shaded areas 
share the same symmetries with the triangular magnetic 
Brillouin zone. We have noted that as J/t — > they 
shrink to the magnetic Goldstones modes. Conversely, 
when J/t > 2.5 there is a finite QP weight anywhere and 
therefore, the full BZ becomes shaded. For the particu- 
lar value J/\t\ = 2, the QP weight vanishes only around 
the M point, as it is displayed by the solid line of Fig. 
QJ This figure also shows the strong k dependence of the 
QP weight for both t signs. Again, a signature of the 
remnant particle-hole symmetry mentioned above is ev- 
idenced: the dashed line (t < 0) resembles the solid one 
(t > 0) shifted by the wave vector Q. In a previous work 
the hole dynamics in a triangular antiferromagnet has al- 
ready been studied using the SCBA^i, but the vanishing 
of the quasiparticlc weight for positive t has been over- 
looked. Our findings, however, are supported by a de- 
tailed and careful treatment of the parameters involved: 
the convergence factor ry, the extrapolation to the ther- 
modynamic limit of the cluster sizes N, and the number 
of frequencies ujliL21. 

By computing the quasiparticle energy scaling with J/t 
it is possible to get some insight of the magnetic cloud 
surrounding the holoii It is known that for a ferromag- 
netic polaron the QP energy scales as e ~ (J/t) - 5 , while 
for an AF polaron it follows a e ~ (J/t) ' 66 dependence. 
In our case we have estimated a e ~ (J/t) a scaling, being 
a ~ 0.55 (~ 0.65) for positive (negative) t in the range 
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FIG. 4: Quasiparticle weight along the T — M' — K — M path 
(see Fig. |HJ Solid line is for positive t and dashed line for 
negative t . 
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FIG. 5: Spectral function for the F, M',K, M points of the 
Brilluoin zone for J/\t\ = 0.4 and negative t. 



J/\t\ < 0.1. A possible interpretation of such results 
is the following: the non-collinear short range magnetic 
correlations around the hole can be thought as a super- 
position of ferro and AF correlations, and then for t < 
the local environment around the hole enhances its AF 
character, while for t > the ferromagnetic character is 
favored. It should be noticed that the opposite behavior 
holds in a simple toy model consisting of a three-sites 
t — J model with one holeiS. Actually, we have checked 
that for cluster sizes smaller than TV = 21, this behav- 
ior is reversed. For instance, for N = 12, a ~ 0.57 for 
positive t while for negative t it turns out a ~ 0.44. 



5 



C. Spectral function structure 

In this section we discuss the main features of the spec- 
tral functions in the thermodynamic limit. The arising 
structures can be traced back to the two hole motion 
processes present in the effective Hamiltonian @. The 
spectral functions at the T, M', K and M points are dis- 
played in Fig. for J/\t\ = 0.4 and negative t . The 
spectra extend over a frequency range of order 9i; that 
is, the non-interacting electronic energy scale. It is pos- 
sible to discern two distinct regions: a low energy sector 
of magnetic origin (magnon-assisted) , and a high energy 
sector with a broad resonance related to the free hop- 
ping process. This resonance, located at an energy near 
ek, has a finite lifetime of order ~ 4 J and a dispersion of 
order t . Furthermore, the low energy sector consists of a 
quasiparticle peak and higher energy resonances with a 
narrow dispersion of order J as the momentun varies. We 
will show below that the latter resonances can be identi- 
fied with string excitations. At V and M' both sectors are 
well separated in energy, while at K and M they overlap. 
From r to M there is a spectral weight transfer from the 
high energy sector to the QP peak as it can be also no- 
ticed quantitatively in the lower panel of Fig. [3 In Fig. 
El we show the magnetic resonances above the QP peak 
for J/\t\ =0.1 (at the ground state momentum M). We 
choose this particular parameter value in order to make 
them more visible, although they are clearly noticeable 
for J/t < 0.5 . In the inset it is shown the energies 
of the first three resonances versus (J/t) 2 / 3 . The linear 
dependence of their energies allow us to associate these 
resonances with string excitations of a hole confined by 
an effective linear potential of order As we have 

previously stated, there is a dynamical enhancement of 
the AF environment around the hole when t is negative, 
and this effect is the responsible for the strong confining 
potential seen by the hole, like in the square lattice case. 

For positive t (Fig. 0, the low energy sector consists 
of a QP peak, when it exists (r and K), and there is no 
signature of string-like excitations. This result is a man- 
ifestation of the ferromagnetic-like environment around 
the hole, pointed out before, for this t sign. Again the 
high energy sector has a broad bandwith of order ~ 9t. 
It is important to notice that at the ground state mo- 
mentum r the QP is very robust, being the only relevant 
feature of the spectral function. Otherwise, at M 1 and 
M the spectral functions are almost structureless. 



IV. CONCLUDING REMARKS 

We have studied the effect of geometrical frustration 
on a single hole moving in an antifcrromagnetic back- 
ground. We chose the triangular t — J model which 
at half filling presents a non-collincar magnetic ground 
state. Using the spinless fermion representation and the 
self-consistent Born approximation we derived an effec- 
tive Hamiltonian with two different hole motion mech- 
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FIG. 6: Spectral functions at the M point for J/\t\ = 0.1 
and negative t. Inset: Energy scaling versus (J/t) 2 ' 3 of the 
first three peaks labeled as I, II and III. 




FIG. 7: Spectral function for the T,M',K,M points of the 
Brilluoin zone for J/\t\ = 0.4 and positive t. 



anisms. One of them arises due to spin-flip processes 
(magnon-assisted) and the other one has its origin from 
the ferromagnetic components of the non-collinear mag- 
netic order. The latter allows the free hole motion with- 
out emission or absorption of magnons and it is absent in 
collinear antiferromagnets. The reliability of the analyt- 
ical method (SCBA) is supported by the excellent agree- 
ment with the exact diagonalization method on small size 
clusters. 

Unlike previous worksi&2I, we have performed a de- 
tailed description of the whole structures of the spectral 
function and we were able to trace them back to the two 
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hole motion processes mentioned above. In particular, 
we have done a careful analysis of quasiparticles, strings, 
and free hole resonances in the thermodynamic limit. 

We have found new and reliable results of the low en- 
ergy sector of the spectra. Our main result is the vanish- 
ing of the quasiparticle excitation for a positive integral 
transfer t at some momenta and at finite J. Also, we have 
uncovered other important change of the spectra under t 
sign reversal: the presence of string-like excitations only 
for t negative. We suggest that these features are a con- 
sequence of different magnetic environment around the 
hole tuned by the t sign. Based on the quasiparticle en- 
ergy scaling with J/i, we argue that for negative t the 
local environment around the hole enhances its AF char- 
acter, while for the other sign the ferromagnetic charac- 
ter is favored. For small clusters (N < 12), this behavior 
is reversed, in agreement with the three-sites problem 
worked out by Koretsune and Ogata^S A deeper study 
of the magnetic environment using the hole wave function 
within the SCBA will be presented elsewhere^. 

The dynamics of a hole in the triangular lattice, 



for negative t, exhibits the same features of the well 
known square lattice cas o 12 ' 19 : a well defined quasiparti- 
cle, string excitations, and antiferromagnetic correlations 
around the hole. On the contrary, for positive t the dy- 
namics is more irregular: no well defined quasiparticle, no 
string excitations, and ferromagnetic correlations around 
the hole. 

For decades researchers have tried to look for spin- 
liquid phases that upon hole doping, give rise to non- 
conventional excitations 1 . Now we give firm evidence 
that such non-conventional excitations can be found in 
non-collinear spin-crystal phases like the one present in 
the triangular antiferromagnet. 
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